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The self-consistent Green's function method is applied to '^ O using a G-niatrix and Vucom as 
effective interactions, both derived from the Argonne uig potential. The present calculations are 
performed in a larger model space than previously possible. The experimental single particle spectra 
obtained with the G-matrix are essentially independent of the oscillator length of the basis. The 
results shows that Vucom better reproduces spin-orbit splittings but tends to overestimate the gap 
at the Fermi energy. 
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A fundamental problem in nuclear physics is how to 
obtain descriptions of finite nuclei starting from a mi- 
croscopic nuclear Hamiltonian. Much progress has been 
achieved for few body systems. The Green's function 
Monte Carlo jj technique is able to give exact results 
up to A=12, while the no-core shell model [3] is accurate 
for even larger nuclei. A wide range of exact methods 
is also available for very light systems Q- Other recent 
attempts to push the limits of ah initio methods into 
the medium mass region have focussed on the nucleus of 
^^O and its neighbor isotopes ;4i,,5,]. These gave predic- 
tions for separation energies and spin orbit splitting of 
the orbits near the Fermi level. Coupled cluster theory 
appears to produce converged results for these nuclei Q ■ 
These achievements have been possible by computing the 
contributions of long-range correlations (LRC) directly 
within very large models spaces where, however, one still 
needs to employ a proper effective interaction that ac- 
counts for the excluded degrees of freedom. In particu- 
lar the effects due to short-range correlations (SRC) can 
be separated efficiently by such partitioning procedure, 
since they are characterized by high momenta degrees 



freedom 
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In ah initio methods the effective interaction is to 
be computed microscopically. Typically, two classes 
of approaches are possible to derive effective interac- 
tions from a standard realistic nucleon-nucleon force. 
Block-Horowitz theory makes use of the Feshbach pro- 
jection formalism to devise an energy dependent inter- 
action 0, Q. This gives solutions for every eigenstate 
with non zero projection onto the model space, however, 
the energy dependence severely cor apli cates the calcu- 
lations. The G-matrix interaction [l^ is a first order 
approximation to the Block-Horowitz scheme. Alterna- 
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tively, one can employ a proper unitary transformation 
to map a finite set of solutions of the initial Hamilto- 
nian into states belonging to the model space. In this 
case, one has the advantage to work with an energy in- 
dependent interaction. Examples of such approach are 
the Lee-Suzuki method 11], Yiow-k [l3 and the uni- 
tary correlator operator method (UCOM) 13, 1% \W{ . 
The UCOM formalism is such that one can apply the 
inverse transformation to reinsert SRC into the nuclear 
wave function. A discussion of the similarities and differ- 
ences between Lee-Suzuki and Block- Horowitz is given in 
Ref. [lal- These methods, in principle, generate effective 
many-nucleon forces in addition to the two-nucleon (2N) 
interaction. In practice, however, calculating medium 
and large nuclei one is forced to work with only 2N Hamil- 
tonians (or at most with weak 3N forces). It is therefore 
important to investigate how this truncation affects the 
results using the different approaches outlined above. 

In Ref. |T3| we proposed to employ a set of Fad- 
deev equations within the self-consistent Green's func- 
tion (SCGF) approach Q to obtain a microscopical de- 
scription of LRC. This allows to couple simultaneously 
quasiparticles to both particle-hole (ph) and particle- 
particle/hole- hole (pp/hh) collective excitations. Such 
formalism was later applied to ^^O to investigate mech- 
anisms that could possibly quench the spectroscopic fac- 
tors of mean field orbits jl^l- These calculations were 
already performed in a no-core fashion. However, the 
model space employed was still somewhat limited and 
phenomcnological corrections were applied to tune the 
values of specific single particle (sp) energies (that allows 
to study correlations by artificially suppressing the cou- 
plings among selected excitation modes). Note that here 
and in the following we use the terms sp energies and 
sp spectra to refer to the poles of the one-body Green's 
function (defined below Eq. ^ ). These represent the ex- 
citation energies of the A±l neighbor nuclei, which are 
observable quantities. In this letter the calculations of 
Ref. [l8J are repeated by avoiding any phenomenology 
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FIG. 1: Example of a Feynman diagram included in the all- 
order summation generated by the set of Faddeev equations. 
Double lines represent the dressed one-particle Green's func- 
tion g{Lo), which propagates quasiparticles (rightward arrows) 
and quasiholes (leftward arrows). The ellipses propagate col- 
lective excitations of the nucleus (Eqs. l(HJ and @ ). 



and employing a large model space. We discuss the re- 
sults of 2N interactions belonging the two types discussed 
above, namely a standard G-matrix and Vucom- 

We consider the calculation of the sp Green's function 
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from which both the one-hole and one-particle spectral 
functions, for the removal and addition of a nucleon, can 
be extracted. In Eq. (UJ, A:"^ = (^^+^141^^) {y^ = 
(^fe ^I'^al^o^)) ^'^^ ^^^ spectroscopic amplitudes for the 



excited states of a system with A + 1 (^ — 1) particles 



and the poles e+ = ^n^^ 
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correspond to the excitation energies with respect to the 
v4-body ground state. The one-body Green's function 
can be computed by solving the Dyson equation [I^ l2y| . 
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where the irreducible self-energy Yj*g{u}) acts as an effec- 
tive, energy-dependent, potential that governs the single 
particle behavior of the system. The self-energy is ex- 
panded in a Faddeev series as in Fig. Q] which couples 
the exact propagator gap{uj) (which is itself a solution 
of Eq. (0) ) to other phonons in the system [l7|. The 
relevant information regarding ph and pp/hh collective 
excitations is included in the polarization and the two- 
particle propagators. Respectively, 
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which describe the one-body response and the propaga- 
tion of two-particles/two-holes. In this work, n(a;) and 
g^^{io) are obtained by solving the dressed RPA (DRPA) 
equations [2ll |23 , which account for the redistribution 
of strength in the sp spectral function. Since this infor- 
mation is carried by the correlated propagator gap{Lo)^ 
Eq. (O, the SCGF formalism requires an iterative solu- 
tion. It can be proven that full self-consistency guaran- 
tees to satisfy the conservation of the number of particles 
and other basic quantities |23l| . 

The coupled cluster studies of Ref. |a, |2J| found that 
eight major harmonic oscillator shells can be sufficient 
to obtain converging results for ^^O with G-matrix in- 
teractions. At the same time, the experience with the 
calculations of Ref. [131 suggests that high partial waves 
do not contribute sensibly. In this work, all the orbits 
of the first eight shells with orbital angular momentum 
I < 4 were included. Inside this model space a G-matrix 
and the Vucom potential were employed as effective in- 
teractions. The former was computed using the CENS 
library routines [Ifl l2al. For the latter, the UCOM 
matrix-elements code |2(tJ was employed with the con- 
straint /^ = 0.09 fm^. This choice of the UCOM cor- 
relator reproduces, in perturbation theory, the binding 
E"^^^) energies of several nuclei up to ^°*Pb 2J]- In both cases 



the Argonne wig potential |23 was used as starting inter- 
action. However, we chose to neglect the Coulomb and 
the other charge independence breaking terms (i.e., to 
set a — 0) in the present work. In applying the SCGF 
formalism, the Hartree-Fock (HE) equations (Brueckner- 
Hartree-Fock (BHF) for the G-matrix) where ffist solved 
for the unperturbed propagator g)^J (w). This was 
employed in the first iteration. After that, the (dressed) 
solution gap{^) was used to solve the DRPA for n(a;) and 
g^^ijj) and then the Faddeev equations for an improved 
self-energy. In each iteration the two most important 
fragments close to the Fermi level were retained for each 
partial wave, both in the quasiparticle and the quasihole 
domains. The remaining strength do not affect sensibly 
the sp states that will be discussed below and was col- 
lected, at each iteration, according to the corresponding 
(B)HF orbitals. In the present work, calculations were it- 
erated until reaching convergence (to within 200 keV) for 
the sp energies nearby the Fermi level. With the available 
computational resources, self-consistence was achievable 
for these states except in the case of the G-matrix with 
harmonic oscillator length hno < 1.8 fm (as discussed be- 
low) . To test the iteration procedure we computed the to- 
tal number of particles obtained with Vucom ■ The first 
calculation, based on the HE propagator, gives A=16.4 



(a 2.5% error). At self-consistency 15.99<A<16.02 (due 
to numerical errors), showing the adequacy of our ap- 
proach. More details on the SCGF/Faddeev formalism 
are given in Ref. 0, IitL lia |. As already noted, how- 
ever, no phcnomenological corrections were applied in the 
present work. 
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FIG. 2: Single particle spectrum obtained with the G-matrix 
as a function of the oscillator length. The dashed lines refer 
to values of buo for which self-consistence was not reached. 

Calculations have been performed for oscillator lengths 
in the interval 6//o = l-6-2.1 fm. The sp spectra obtained 
at self-consistency are shown in Figs. El and El In the 
case of a G-matrix, these orbitals appear independent, 
within numerical accuracy, of the oscillator frequency for 
bi/o > 1.9 fm. Below this the trend is similar. How- 
ever, the structure of the quasiparticle spectral function 
becomes increasingly complex for the l—l waves that cor- 
respond to the pf shell. These break in more than two 
relevant fragments and their sp energies depend strongly 
on bno, suggesting a non negligible contamination of the 
c.o.m. motion. In this situation the above approach for 
obtaining self-consistency cannot be carried out reliably. 
We show the result of a typical iteration by dashed lines 
in Fig. El but discard these results in the following dis- 
cussion. This complication does not happen for Vjjcom, 
which is a softer interaction and generates, for each shell, 
at most one main fragment and a smaller satellite peak 
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FIG. 3: Self-consistent single particle spectrum obtained with 
VucoM as a function of the oscillator length. 
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Vc/com: 

A£p,/,_P3/, 4.7 4.4 4.1 4.5 4.4 4.1 

AEd.-^^^^d,,^ 4.9 4.4 4.1 3.9 3.6 3.3 



TABLE I: Spin orbit splittings (in MeV). 



near the Fermi energy. Fig. El shows that the spin orbit 
splittings for this interaction are approximately constant, 
although the sp energies are not yet independent of the 
oscillator length. This can be understood considering 
that these spin orbit partners correspond to particularly 
simple and similar configurations (one particle or one hole 
on top of the correlated ground state) . Conversely, sep- 
aration energies are linked to the to total hinging energy 
of neighbor isotopes. Larger model spaces will probably 
be required for a full convergence with Vucom ■ 

The splittings obtained from both interactions are re- 
ported in Tab. ^ For the Op orbits these are practically 
constant for all the oscillator lengths. In general the G- 
matrix predicts about a half of the experimental value. 
Better solutions are obtained with the present choice of 
the UCOM correlator. For the Od orbits the results for 
the two interactions are more similar to each other but 
not totally independent of the oscillator length. 

The same type of LRC studied here (which are pre- 
dominantly of 2plh and 2hlp in character) were also 
considered in ReL p|. There, the effective interaction 
was derived in the unitary-model-operator approach and 
an explicit diagonalization was performed. The result- 
ing sp energies and spin-orbit splitting, however, showed 
a stronger dependence on the oscillator length than the 
one found in this work. The spectrum obtained in Fig. El 
is nearly convergent, suggesting that the all order sum- 
mation employed here and the proper accounting of the 
fragmentation of the sp strength allow to select relevant 
configurations beyond the bare 2plh/2hlp level. Cou- 
pled cluster calculations are also available for ^^O with an 
analogous wig/G-matrix "q. These authors report split- 
tings of the Op and Od orbits larger than those of Tab.Qby 
about 1.5 and 0.5 MeV, respectively. The present work 
includes LRC only in the form of coupling to small am- 
plitude excitations of the core — which can be described 
at the DRPA level. More complex collective modes are 
also present [231 and should be included for a full solution 
of the many-body problem. We note, for example, that 
the phcnomenological studies of Ref. |l^ suggest further 
contributions to the p^i2 quasihole wave function coming 
from couplings to the first excited O"*" state in oxygen. 
Testing this conjecture would first require being able to 
reproduce the correct excitation energy of this level — 
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FIG. 4: Single particle spectra generated by the G-matrix in 
the HF approximation, after the first iteration and in the full 
self-consistent solution. The points on the left represent the 
experimental values for the 017/015 case. 





lOr 




5r 


> 

1 


Or 

-5r 


d. 


-lOr 


tq 


-15r 




-20 r 




-25 - 



HF 



1 itr. 



All 



,-1/2 



SCGF 



FIG. 5: Same as Fig. 0|but for the Vucom interaction. 



since it can couple effectively only when it is low enough 
in energy. To our knowledge this is still a challenge for 
the available ab initio methods. 

Figures0]and[Slshow the effects of LRC on the sp spec- 
trum for bno— 1-9 fni, and compare to the experimental 
values for the addition/removal of a neutron. For both 
interactions the coupling to collective phonons reduces 
the splitting of the Op orbits, with respect to the HF ap- 
proximation. Including the effects of fragmentation tends 
instead to compress the sd shell and to lower the whole 
spectrum. The self-consistent results for the energy gap 
between particle and hole states, AEp = ^dg ,2 ~Spi/2 ' ^^^ 
13.0 MeV with the G-matrix and 15.4 MeV with Vucom- 
Both of them exceed the experimental value of 11.5 MeV. 
The mean square radii obtained are r„„.s=2.63 fm (G- 



matrix) and r,.„s=2.45 fm (Vucom)- 

Both the UCOM and the renormalization group ap- 
proaches have the capability of shifting the binding ener- 
gies for A=3,4 systems along the Tjon line i30|,|31|. In the 
first case this can be achieved by varying the correlator 
in the tensor-isoscalar channel. For Viow-k, one modifies 
the momentum cut off. Usually, tuning the binding ener- 
gies to the experimental value increases the non locality 
of the interaction and leads improved spin orbit split- 
tings, as seen in Tab.Q] On the other hand, our Vucom 
result for AEp -with 2N forces- overestimates the ex- 
periment, more than the G-matrix does. This behavior 
is seen already at the HF level for soft interactions like 
Vucom and Viow-k |2^il^ and it is only slightly modi- 
fied by the LRC considered here. In all cases, three-body 
forces appear necessary in order to reproduce the whole 
spectrum of observations. We note, however, that the 
UCOM method offers some advantages to reduce the con- 
tribution needed from many-body forces since it allows 
to treat SRC in different channels separately [l^. 

In conclusion, SCGF calculations have been performed 
for the first time in a large model space, including up 
to eight oscillator shells. Long-range correlations in the 
form of coupling sp to ph and pp/hh RPA modes were 
investigated for ^^O. A comparison was made between 
the results of a G-matrix and the Vucom effective inter- 
actions, both derived from same realistic potential (Ar- 
gonne t^is). The spectra of adjacent nuclei were found 
to be nearly convergent for the G-matrix, while they de- 
pend only weakly on the oscillator length for Vucom - 
In general it was found that the LRC effects considered 
here, tend to compress the spectra of A±l nuclei but do 
not affect sensibly the the gap between quasiparticle and 
quasihole energies at the Fermi level. 
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